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eigenvalues of this many-body system and perform the summation over the entire spectrum of 
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determinant with the Airy kernel yielding the universal Tracy- Widom distribution, which is known to 
describe the statistical properties of the Gaussian unitary ensemble as well as many other statistical 
systems. 
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I. INTRODUCTION 

Directed polymers in a quenched random potential have been the subject of intense investigations during the 
past two decades (see e.g. [![). Diverse physical systems such as domain walls in magnetic films Q, vortices in 
superconductors wetting fronts on planar systems Q, or Burgers turbulence [f| can be mapped to this model, 
which exhibits numerous non-trivial features deriving from the interplay between elasticity and disorder. In the 
most simple one-dimensional case we deal with an elastic string directed along the r-axis within an interval [0,L]. 
Randomness enters the problem through a disorder potential V[</>(t),t], which competes against the elastic energy. 
The problem is defined by the Hamiltonian 

H^(t), V] = £ dr{ ~ [dMr)] 2 + V[4>(t), t] } ; (1) 



where in the simplest case the disorder potential V[<f>, r] is Gaussian distributed with a zero mean V(4>,t) = and 
the ^-correlations: 



V{4>, t)V{4>', t') = uS{t - t')5(4> - 4>') 



(2) 
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Here the parameter u describes the strength of the disorder. Historically, the problem of central interest was the 
scaling behavior of the polymer mean squared displacement which in the thermodynamic limit {L —y oo) is believed 
to have a universal scaling form (cj) 2 )(L) oc L 2< » (where (...) and (...) denote the thermal and the disorder averages), 
with £ = 2/3, the so-called wandering exponent. More general problem for all directed polymer systems of the 
type, Eq.((T]), is the statistical properties of their free energy fluctuations. Besides the usual extensive (linear in L) 
self-averaging part foL (where /o is the linear free energy density), the total free energy F of such systems contains 
disorder dependent fluctuating contribution F, which is characterized by non-trivial scaling in L. It is generally 
believed that in the limit of large L the typical value of the free energy fluctuations scales with L as F oc L 1 / 3 (see 
e.g. [6-9]). In other words, in the limit of large L the total (random) free energy of the system can be represented as 

F = f Q L + cL 1 ' 3 f (3) 

where c is the parameter, which depends on the temperature and the strength of disorder, and / is the random 
quantity which in the thermodynamic limit L — > oo is described by a non-trivial universal distribution function P*(/). 
The derivation of this function for the system with 5-correlated random potential, Eqs.([T])-([2]) is the central issue of 
the present work. 

For the string with the zero boundary conditions at r = and at r = L the partition function of a given sample is 

Z[V] = / T>[<t>(T)] e~^' y] (4) 

where f3 denotes the inverse temperature. On the other hand, the partition function is related to the total free energy 
F[V] via 

Z[V] = exp(-/?F[F]) (5) 

The free energy F[V] is defined for a specific realization of the random potential V and thus represent a random 
variable. Taking the iV-th power of both sides of this relation and performing the averaging over the random potential 
V we obtain 



Z N [V] = Z[N,L] = exp{-(3NF[V}) (6) 

where the quantity in the lhs of the above equation is called the replica partition function. Substituting here F — 
faL + cL 1 / 3 f, and redefining 

Z[N, L] = Z[N, L] e~ pNfoL (7) 

we get 

Z[N,L]=exp(-\Nf) (8) 

where A = fid 1 / 3 . The averaging in the rhs of the above equation can be represented in terms of the distribution 
function Pl(I) (which depends on the system size L). In this way we arrive to the following general relation between 
the replica partition function Z[N,L] and the distribution function of the free energy fluctuations Pl(/): 

/+oo 
dfp L (f)e- XN f (9) 
-oo 

Of course, the most interesting object is the thermodynamic limit distribution function P*(f) — lim£_ i . 00 Pl{/) which 
is expected to be the universal quantity. The above equation is the bilateral Laplace transform of the function Pi(/), 
and at least formally it allows to restore this function via inverse Laplace transform of the replica partition function 
Z[N, L]. In order to do so one has to compute Z[N,L\ for an arbitrary integer N and then perform analytical 
continuation of this function from integer to arbitrary complex values of N. In Kardar's original solution [9j, after 
mapping the replicated problem to interacting quantum bosons, one arrives at the replica partition function for 
positive integer parameters N > 1. Assuming a large L — > oo limit, one is tempted to approximate the result by 
the ground state contribution only, as for any N > 1 the contributions of excited states are exponentially small for 
L — > oo. However, in the analytic continuation for arbitrary complex TV the contributions which are exponentially 
small at positive integer N > 1 can become essential in the region N — > 0, which defines the function P(f) (in other 
word, the problem is that the two limits L — > oo and N — > do not commute da mi). As a consequence, the resulting 
distribution function P*(f) exhibits severe pathologies such as the vanishing of its second moment, which assumes that 
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the distribution function is not positively denned. Nevertheless, in terms of this approximation (assuming he universal 
scaling L 1 / 3 of the free energy fluctuations) one can derive the left tail of the distribution function P*(/) which is 

given by the asymptotics of the Airy function, P*(/ — > — oo) ~ exp(— 1|/| 3 ^ 2 ) [ljj]. For the first time the form of the 
right tail of this distribution function, — > +00) ~ exp[— (const) f 3 ^ , has been derived in terms of the optimal 
fluctuation approach (l3| . where it has been demonstrated that both asymptotics (left and right) of the function 
P*(f) are consistent with the Tracy- Widom (TW) distribution Originally the TW distribution has been derived 
in the context of the statistical properties of the Gaussian Unitary Ensemble while at present it is well established 
to describe the statistics of fluctuations in various random systems [l5l - [2Cl | , which were widely believed to belong to 
the same universality class as the present model [TH, Hi], [22| ■ Finally, quite recently, the one-point distribution of the 
solutions of the KPZ-equation (which is equivalent to the present model) has been derived for an arbitrary value of 
L, and it has been shown that in the limit L — > 00 this distribution turns into the Tracy- Widom distribution [23l |24|. 

In the recent paper [25[ an attempt has been made to derive the free energy distribution function P*(/) via the 
calculation of the replica partition function Z[N,L\ in terms of the Bethe-Ansatz solution for quantum bosons with 
attractive ^-interactions which involved the summation over the entire spectrum of exited states. Unfortunately, the 
attempt has failed because on one hand, the calculations contained a kind of a hidden " uncontrolled approximation" , 
and on the other hand, the analytic continuation of obtained Z(N, L) was found to be ambiguous. Later on it turned 
out that it is possible to bypass the problem of the analytic continuation if instead of the distribution function itself 
one would study its integral representation 



W(x) 



df P*(f) 



(10) 



which gives the probability to find the fluctuation / bigger that a given value x. Formally the function W(x) can be 
defined as follows: 



W(x) 



lim 



E 

N=0 



(-1) 

Nl 



N 



■exp(\Nx)Z N = lim exp[-exp(A(x - /))] = 6{f - x) 



(11) 



On the other hand, the probability function W(x) can be computed in terms of the replica partition function Z[N, L] 
by summing over all replica integers 



W{x) 



lim > 



-1 



JV=0 



AT! 



• exp(XNx) Z[N, L] 



(12) 



In terms of this approach, quite recently it has been reported that performing the summation over the entire spec- 
trum of excited states of TV- particle quantum boson system the above series, eq. (fT2"|) . can reduced to the Fredholm 
determinant with the Airy kernel which is known to yield the Tracy- Widom distribution [2H |27| . 

Present paper provides the detailed derivation of this result. In particular, Appendices A and B contain detailed 
analysis of the properties of ./V-particle wave functions of quantum bosons with both repulsive and attractive 5- 
interactions. Besides, basing on the definition of the Fredholm determinant with the Airy kernel, in Appendix C it is 
demonstrated how the explicit form of the Tracy- Widom distribution can be derived in terms of the solution of the 
Panlcvc II differential equation. 

II. MAPPING TO QUANTUM BOSONS 

Explicitly, the replica partition function, Eq.©, of the system described by the Hamiltonian, Eq.(fT|), is 



N „0 a (L) =O 

Z{N,L) = T[ / P0 Q (r)exp 
a=1 J«Mo)=o 



P I drJ2{l[drMr)] 2 + V[cP a (r),T]} 



(13) 



Since the random potential V[<fi, r] has the Gaussian distribution the disorder averaging (...) in the above equation is 
very simple: 



exp 



pL N - I" 2 L N 

-13 / dTj2v[Mr),r] = exp tL / / drdr' V F[^(r), T]V\<k{r>), t'\ 

Ja 0=1 J L J J ° 0,6=1 



(14) 
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Using Eq.© we find: 

Z(N,L) = TT / Xty (r) exp --j3 / drl V[c> t ^(t)] - ^ V *[*„(t) - &(r)] } 



a, 6=1 



(15) 



It should be noted that the second term in the exponential of the above equation contain formally divergent contribu- 
tions proportional to 5(0) (due to the terms with a = b). In fact, this is just an indication that the continuous model, 
Eqs.|T])-([2]) is ill defined as short distances and requires proper lattice regularization. Of course, the corresponding 
lattice model would contain no divergences, and the terms with a = b in the exponential of the corresponding replica 
partition function would produce irrelevant constant ^L/3 2 uNS(0) (where the lattice version of (5(0) has a finite value). 
Since the lattice regularization has no impact on the continuous long distance properties of the considered system this 
term will be just omitted in our further study. 

Introducing the A^-component scalar field replica Hamiltonian 



1 f L ( N N 
H N [ct>] = ^ dri^idrMr)] 2 -pu^SiM^-Mr)]) (16) 

*' Q \o=l a^b 



for the replica partition function, Eq.([T5|). we obtain the standard expression 

N „0 a (L)=O 

Z(N,L) = T\ V<j> a { T ) e-f" 1 *™ (17) 

a=1 i0 a (O)=O 

where <j> = {<f>i, . . . , 4>n}- According to the above definition this partition function describe the statistics of N 
(^-interacting (attracting) trajectories <j> a (T) all starting (at r = 0) and ending (at r = L) at zero: </> a (0) = (f> a (L) = 
In order to map the problem to one-dimensional quantum bosons, let us introduce more general object 

*(x;t) = TT/ -DM?) e- pHlTl * ] (18) 

a=l J M0)=0 

which describes N trajectories <j)a(j) all starting at zero (<^ a (0) = 0), but ending at r = t in arbitrary given points 
{xi, ...,xjy}. One can easily show that instead of using the path integral, ^(x;t) may be obtained as the solution of 
the linear differential equation 

N N 

ft*(x;t) = — + -(3 2 uJ2S(xa-x b )y(x;t) (19) 

" 0=1 a=£b 

with the initial condition 

\£-(x;0)=nf =1 £(z a ) (20) 
One can easily see that Eq. (|19p is the imaginary-time Schrodinger equation 

-at#(x;*) = .#*(x;t) (21) 

with the Hamiltonian 

N N 

V=-^/Z d t- <^ - x b ) (22) 

P a=l a=jtb 

which describes N bose-particles of mass /3 interacting via the attractive two-body potential —j5 2 u5(x). The original 
replica partition function, Eq. (|17|) . then is obtained via a particular choice of the final-point coordinates, 

Z(JV,L) = *(0;L). (23) 

The eigenfunction of A^-particle system of repulsive bosons (u < 0), Eq. (|22p . have been derived by Lieb and Liniger 
in 1963 |28j (for details see Appendix A, as well as Refs. [29, 30]). The spectrum and some properties of the 
eigenfunctions for attractive (u > 0) one-dimensional quantum boson system have been derived by McGuire [3l| and 
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by Yang [32[ (see also Ref. [33j, |34|). Detailed structure and the properties of these wave functions are described in 
Appendix B. A generic eigenstate of such system consists of M (1 < M < N) "clusters" {ft a } (a — 1, M) of bound 
particles. Each cluster is characterized by the momentum q a of its center of mass motion, and by the number n a 
of particles contained in it (such that X)a=i Ua = ^0' Correspondingly, the eigenfunction VPq^n {x\, ■■■jXn) of such 
state is characterized by M continuous parameters q = (qi, qm) and M integer parameters n = (n\, ...,Um) (see 
Appendix B2, Eq. (jB.28jl ). The energy spectrum of this state is 

1 M K 2 M 

E M {<\,n) = n a ql- — ^(nl-n a ) (24) 

" a=l " a=l 

where 

k = f3 3 u (25) 

A general time dependent solution \I/(x, t) of the Schrodinger equation (fll?]) with the initial conditions, Eq. ([2H|) . can 
be represented in the form of the linear combination of the eigenfunctions ^q^n (x) (Appendix B4, Eq. (IB.5ip ): 

N / , 

*(x,t) = J2 E / ^ *^(*)*^>) exp[-£ M (q,n)i] (26) 

M=l n J 

where the summations over the integer parameters n a and the integrations over the momenta q a are performed in a 
restricted subspace, Eqs. (|B.43l) - (|B.46[) and (|B.52I) . which reflects the specific symmetry properties of the eigenfunctions 

^qjn' (x). Correspondingly, according to Eq. (|23|) . for the replica partition function of the original directed polymer 
problem one gets (Appendix B4, Eq. (|B.60l) ') 

N 1 

Z(N,L) = y — 

Af=l 

where due to the symmetry of the function /(q, n) = j^q^n (0)| Z exp[— Em(<1, n) L] with respect to permutations of 
all its M pairs of arguments {q a ,n a ) the integrations over M momenta q a can be extended to the whole space Rm 
while the summations over n a 's are bounded by the only constrain n a = N (for simplicity, due to the presence 

of the Kronecker symbol S(J2 a n a , Nj, the summations over n Q 's are extended to infinity). 



n 



i-oo 



dq a 
2tt 



E 



^2 n a , 



n) |* q ^(o)| 2 e- 



-E M (q,n)L 



(27) 



III. REPLICAS PARTITION FUNCTION AND THE FREE ENERGY DISTRIBUTION FUNCTION 



Using the explicit form of the wave functions "fq.n (x), Eq. (|B.28[) . the expression in Eq. (|27|) for the replica partition 
function can be reduced to (Appendix B4, Eq. (|B.62|) - (jB.63|l ) 



z(n,l) =e 



Z(N,X) 



(28) 



where / = ^(3 4 u 2 — ln(/3 3 w) is the linear (selfaveraging) free energy density (cf. Eq.©, and 



Z(N.L) = N\ 



dq 



1 



2ttkN 

M oo 



exp 



2(3 q 24/3 



M=2 ' L a=l n a =l 



dq a 



M 



2iTKn c 



M 



M 



a<f) 



M 



|?« - qp - ^(na - np)\ 
\q a - qp - % -r{n a + np)\ 



x exp 



L ^ — > q K 2 L ^ — > 

E n ^ a + ^l^ n a 



20 



24(3 ^ 

' a = l 



(29) 



The first term in the above expression is the contribution of the ground state (M — 1), while the next terms (M > 2) 
are the contributions of the rest of the energy spectrum. 
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The terms cubic in n a in the exponential of Eq. (|29f can be linearised with the help of Airy function, using the 
standard relation 



1 f + °° 
exp(-A 3 ?i 3 ) = / dy Ai(y) exp(An) 

3 J-oo 

1 /3 

Redefining the momenta, g Q = [Bk/Ia p a and introducing a new parameter 



A(L) 



1 / i 



2V/3 



1/3 



i(/>Vi) 



1/3 



(30) 



(31) 



after shifting the Airy function parameters of integration y a — > y a +p 2 the expression for Z(N, A) becomes sufficiently 
compact: 



Z{N,X) = Nl 



dydp 
4ttAA 



Ai(y + p 2 )e 



2\ pAiVy 



(32) 



iV r M 



dy a dp a 2\C^ n aVa 

— Ai(y a +p a )e 



M 

n 



|A(n Q - np) - «(p Q 



\X(n a +np) —i(p a —pp)\ 



M 



Now, using the Cauchy double alternant identity 



Ylv<p{ a <x -ap)(b a -bf)) 



n^=iK -bp) 

the product term in eq. (|32l) can be represented in the determinant form: 



= (-lj^-^det 



= 1,...M 



M 

n 



|A(n Q - rip) - i(p a -pp)\ 



a</3 \H n a +np) -i{p a -Pp)\ 



M 



a=l 



det 



Xi2 a — ip a + Xnp + ipp 



a,P=l,...M 



« E n «. N 



(33) 



(34) 



Substituting now the expression for the replica partition function Z{N 1 A) into the definition of the probability function, 
eg. (1121) . we can perform summation over N (which would lift the constraint Yl a =i n a = N) and obtain: 



(-1) 



M 



M=l 



Ml 



M 

n 



-HOC: 



dy a dp a 
2vr 



Ai(y Q +^) £ (-lr^e 



n a - 1 pAn a (i/c +x) 



det 



f 



-An Q — ip a + Xnp + ipp 
(35) 

The above expression in nothing else but the expansion of the Fredholm determinant det(l — K) (see e.g. [35]) with 
the kernel 



K = K[{n,p);(n',p')] 



+ 00 



d 2 /Ai(y+p 2 )(-l)"- 1 e An(H+x) 



f 



Xn — ip + Xn' + ip' 



Using the exponential representation of this determinant we get 

oo ^ 

W(x) — lim exp — > — 

V ' A^oo L Z-~t M 
M=l 



TrK 



M 



(36) 



(37) 



where 



TrK 



M 



M 

n 



dy a dp 
2tt 



Ai(y a + P 2 a ) £ (_i)»--ie An -^^ 



n a = l 
1 



Substituting here 



(Xni - ipi + Xn 2 + ip2){Xn 2 - ip 2 + A?i 3 + ip 3 )...(Xn M - ipM + Ani + ip\) 



1 f 00 

: — — = / duj a exp[-(An Q - ip a + Xn a+1 + ip a+1 )ui a ] 

M - ip a + Xn a+ i + ip a+ i Jo 



(38) 



(39) 
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one can easily perform the summation over n Q 's. Taking into account that 

oo pAz 

lim £(-l)-ie A " z = lim -——^ = 9(z) (40) 

n=l 

and shifting the integration parameters, y a — > y a — x + uj a + uj a -i and uj a — > ui a +x/2, after taking the thermodynamic 
limit, L — > oo (which according to Eq. (|31[) is equivalent to A — > oo) we obtain 



lim Tr if M = TT / dy a / f2L / du , Q Ai(y Q + j£ + Wq + ^ -0 (41) 



where by definition it is assumed that cjq = w^f. Using the Airy function integral representation, and taking into 
account that it satisfies the differential equation, Ai (t) = t Ai(t), one can easily perform the following integrations: 

/•OO n -f" OO 7 />00 

/ dy / ^Ai( 1 /+p 2 +w + w')e* (w - w ' ) = 2- 1 / 3 / rfyAi(2 1 / 3 W + y)Ai(2 1 /V + 2 ;) (42) 

Ai(2 1 /3 W ) Ai'(2 1 /3 W ') _ Ai'(2 1 /3 W ) Ai(2 1 /3 W ') 
u> — Ul' 

Redefining u> a — > o; a 2 -1 / 3 we find 



A- 

wherc 



lim TrK M = / / ... / dujidw 2 —dui M KA{^i, u 2 )K a (w 2 ,uj 3 )...Ka(ojm, ui) (43) 



K , Mju) Ai'(a/) - Ai'(u) Ai(u') 

K A (u,u) = (44) 

UJ — Ul 

is the so called Airy kernel. This proves that in the thermodynamic limit, L — > 00, the probability function W{x), 
cq. pH)) . is defined by the Fredholm determinant, 

W(x) = det[l-K A ] = F 2 {-x/2 2/3 ) (45) 

where Ka is the integral operator on [— x/2 2 / 3 , 00) with the Airy kernel, eq. (|44l) . The function F 2 (s) is the Tracy- 
Widom distribution [l4| 

F 2 (s) = exp(-^ dt(s-t)q 2 {t)^ (46) 

where the function q{t) is the solution of the Panleve II equation, q" = tq + 2q 3 with the boundary condition, 
q(t — > +00) ~ Ai{t) (see Appendix C). 



IV. CONCLUSIONS 



The Tracy- Widom distribution, Eq. (|46[) . was originally derived as the probability distribution of the largest eigen- 
value of a n x n random hermitian matrix in the limit n — > 00. At present there are exists an appreciable list of 
statistical systems (which are not always look similar) in which the fluctuations of the quantities which play the role 
of "energy" are described by the same distribution function F 2 (s). These systems are: the polynuclear growth (PNG) 
model [15|], the longest increasing subsequences (LIS) model [lg, the longest common subsequences (LCS) [lj]], the 
oriented digital boiling model [l8| , the ballistic decomposition model [l9| , the zero-temperature lattice version of the 
directed polymers with an exponential and geometric site-disorder distribution [20| . Now we can add to this list the 
continuous one-dimensional directed polymers with Gaussian (5-correlated random potential. 
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Appendix A 
Quantum bosons with repulsive interactions 



1. Eigenfunctions 

The eigenstates equation for TV-particle system of one-dimensional quantum bosons with ^-interactions is 

N N 



(A.l) 



(where n = /3 3 u). Due to the symmetry of the wave function with respect to permutations of its arguments it is 
sufficient to consider it in the sector 



X\ < X2 < ... < Xn 
as well as at its boundary. Inside this sector the wave function ^(x) satisfy the equation 



1 N 

= -/J25*(x) 



(A.2) 



(A.3) 



which describes N free particles, and its generic solution is the linear combination of N plane waves characterized by 
N momenta {qi, q2, <7tv} = Q- Integrating Eq. (|A.l| over the variable (xj+i — Xi) in a small interval around zero, 
\xi+\ — Xi\ < e — s> 0, and assuming that the other variables {xj} (with j ^ i, i + 1) belong to the sector, Eq. (|A.2j) . one 
easily finds that the wave function ^(x) must satisfy the following boundary conditions: 



= 



(A.4) 



Functions satisfying both Eq. (|A.3[) and the boundary conditions Eq. (IA.4I) can be written in the form 



y qx ... qH (xi,...,x N ) = *W( X ) = c [Y[[ d x a - d Xb + n]) det expfe 



(c,d) = l,...,N 



(A.5) 



where C is the normalization constant to be defined later. First of all, it is evident that being the linear combination 
of the plane waves, the above wave function satisfy Eq. ([A.3p . To demonstrate which way this function satisfy the 
boundary conditions, Eq. (|A.4[) . let us check it, as an example, for the case i = 1. According to Eq. (IA.5l) . the wave 
function ^^(x) can be represented in the form 



*W(x) = -(a X2 -a xl -*)#W(x) 



where 



N 



*^(x) = Cm [d Xl - d Xa + k] [d X2 - d Xa + k])[ Yl [d Xa - d Xb + k] det exp(zg c x d ) 



a— 3 



N 



3<a<b 



(c,d)=l,...,N 



(A.6) 



(A.7) 



One can easily see that this function is antisymmetric with respect to the permutation of x\ and x-i- Substituting 
Ea.(|A.6|l into Eq. (|A.4l) (with i = 1) we get 



(d X2 - d Xl ) 2 - k 2 



fW(x) 



= 



(A., 



Given the antisymmetry of the l.h.s expression with respect to the permutation of x\ and x^ the above condition is 
indeed satisfied at boundary x\ = x%. 
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Since the eigenfunction \I>q '(x) satisfying Eq. (|A.l[) must be symmetric with respect to permutations of its argu- 
ments, the function, Eq. flA.5l) . can be easily continued beyond the sector, Eq. flA.2l) . to the entire space of variables 
{x 1 ,x 2 , ...,x N } € R N , 



N 



*W(x) = Cm [-i(d Xa -d Xb )+iKsgn(x a -x b ) 



det 



exp(iq c x d ) 



(c,d) = l,...,N 



(A.9) 



where, by definition, the differential operators d Xa act only on the exponential terms and not on the sgn(x) functions, 
and for further convenience we have redefined i N ^ N ~ x )/ 2 C — > C. Explicitly the determinant in the above equation is 



det 



exp(iq c x d ) 



{c,d) = l,...,N 



N 



(A.10) 



where the summation goes over the permutations P of N momenta {91,92, ■•-,9at} over N particles {x\,X2, ...,Xn}i 
and [P] denotes the parity of the permutation. In this way the eigenfunction, Eq. (|A.9p . can be represented as follows 

, N s N 

*W(x) = C^(-l)[ p l M][-i(9 Ia -d Xb ) + ik sgn(x a ~x b )\\ exp[iJ2l P aXa] (A.ll) 

P ^a<6 ' a=l 

Taking the derivatives, we obtain 

/ JV s N 

*W(x) = C^(-l) [p] (n[q Pa -q Pb +iKSgn(x a -x b )n exp[z^ q Pa x a j (A.12) 

P V<b ' a=l 



It is evident from these representations that the eigenfunctions \&q (x) are antisymmetric with respect to permuta- 
tions of the momenta q\,...,qN- 

Finally, substituting the expression for the eigenfunctions, Eq. (|A.5p (which is valid in the sector, Eq. (|A.2[) ). into 
Eq. (|A.3p for the energy spectrum we find 



E = 



1 N 



2/3 



(A.13) 



2. Orthonormality 

Now one can easily prove that the above eigenfunctions with different momenta are orthogonal to each other. Let 
us consider two wave functions V&q ^(x) and ty^f (x) where it is assumed that 

qi < q 2 < - < qw (A. 14) 

q[ < q' 2 < ... < q' N 

Using the representation, Eq. (|A.ll[) . for the overlap of these two function we get 

/+00 
^x*W'(x)*W(x) 
-OC 

= |C| 2 ^(-l)[ p W p '] f +OC d N J(f[[ l (d Xa -d Xb )- l Ksgn(x a ~x b )]) expi-if^faxa]) x 

P,P> "'- 00 \<b ' a=l > 

N s N >. 

(A.15) 



JV \ J V X 

H[ _i ( 9 x a -9 Xb ) +iKSgn(x a - x b )] j exp[i'^2q Pa Xa\ > 

a< b ' a=l ' 
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Integrating by parts we obtain 



(A.16) 



/-\-oo 
d N x e X p[-iJ2 q >, Xa ]x 
t-,t- ~°° 0=1 

, N \ N 

Y\_[-i(9 Xa ~ d Xb ) -insgn(x a - x b )] [-i(d Xa - d Xb ) + insgn(x a - x b )] J exp[i^g Pa x a ] 

v a<6 ' o=l 



f.t- "°° a=l V<b ' a=l 

(A.17) 



Taking the derivatives and performing the integrations we find 

N 



/ N \ r +oo n 

*W' (X) *W (X) = \C\^(-V [P]+[P ' ] (m(^-^) 2 + K2 ]) / ^xexp[i^(^- g ;,K 

p.P' V<b ' J -°° 0=1 



p.p' 



| C | 2 J- (_i)[-P]+[P'] TJ [feo _ qpb? + K 2] JJ^)^ - q ' p ,J 



a<b 



N 



(A.18) 



Taking into account the constraint, Eq. (|A.14[) . one can easily note that the only the terms which survive in the above 
summation over the permutations arc P = P', all contributing equal value. Thus, we finally get 



N 



r N 



o=l 



( X )<> (x) = id 2 m n - ^ + « 2 ] ) n & 

With the normalization constant 



(A.19) 



(A.20) 



we conclude that the set of the eigenfunctions, Eq. (|A.ll[) or (|A.12|) . are orthonormal. The proof of completeness 
of this set is given in Ref. [3Cj. It should be noted that the above wave functions present the orthonormal set of 
eigenfunctions of the problem, Eq. ([A.ip , for any sign of the interactions k, e.i. both for the repulsive, n < 0, and for 
the attractive, k > 0, cases. However, only in the case of repulsion this set is complete, while in the case of attractive 
interactions, n > 0, in addition to the solutions, Eq. (|A.lip . which describe the continuous free particles spectrum, 
one finds the whole family of discrete bound eigenstates (which do not exist in the case of repulsion) . 



Appendix B 
Quantum bosons with attractive interactions 



1. Ground state 



The simplest example of the bound eigenstate defined by eq. (|A.1[) (with n > 0) is the one in which all N particles 
are bound into a single " cluster" : 



*W(x) = C exp 



"1 



N I N 

Xa ~ i k 5Z \ xa ~ 



(B.l) 



a=l 



,6=1 
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where C is the normalization constant (to be defined below) and q is the continuous momentum of free center of mass 
motion. Substituting this function in Eq. (|A.ip , one can easily check that this is indeed the eigenfunction with the 
energy spectrum given by the relation 



E 



i N r i N 



(B.2) 



where it is assumed (by definition) that sgn(O) = 0. Since the result of the above summations does not depend on the 
mutual particles positions, for simplicity we can order them according to Eq. (|A.2| ). Then, using well known relations 



N 



6=1 



E sgn{x a - Xb) = -(N+l-2a) 
N 1 



a=l 
N 



for the energy spectrum, Eq. (|B.2j) . we get 

E = 



E fl2 = l(N + l)(2N+l) 



^-^-N)^E l[q ,N) 



The normalization constant C is denned by the orthonormality condition 



*«»*«(x) 



dx x ...dx N ^ } *(x)*W(x) = (27r)S(q-q') 



(B.3) 
(B.4) 
(B.5) 

(B.6) 
(B.7) 



Substituting here Eq. flB.ip we get 

r + oo r N ^ N 

*^ ) *(x)*^ 1) (x) = |Cf / dx x ...dx N exp i(q - q') E x a ~ -k E 1^ 

■/ — OO . 1 * „ l_1 



o=l 



.,6=1 



i-+oo /-+oo 






/ dxi / (2X2... 


. dxN exp 


i(q 


' — oo J a;i 







N 



N 



KQ - <l')J2 x a + kJ2(N + 1 - 2a)x a 



a=l 



(B. 



where for the ordering, Eq. (|A.2[) . we have used the relation 



9 E !*« 



iV 



^6 



E(^V + l-2a)x 



(B.9) 



a, 6=1 a=l 

Integrating first over xjy, then over xjv-i, and proceeding until x\, we find 

,N-\ , x P+ oo 



*J } *(x)vI/«(x) = |C| 2 AT! J] 



/Vic 

|C| 2 ^r(2 7 r)%-,') 



dxi exp [iN(q — q')xi 



N\k n 

According to Eq. (|B.7|) this defines the normalization constant 



< = . c«., (9) 



(B.10) 



(B.ll) 
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Note that the eigenstate described by the considered wave function, Eq. (|B.ip . exists only in the case of attraction, 
k > 0, otherwise this function is divergent at infinity and consequently it is not normalizable. 



It should be noted that the wave function, Eq. (|B.l[) . can also be derived from the general eigenfunctions structure, 
Eq. (IA.12[) . by introducing (discrete) imaginary parts for the momenta q a . We assume again that the position of 
particles are ordered according to Eq. (|A.2p . and define the particles' momenta according to the rule 



(la 



q - -k(N + 1 - 2a) 



(B.12) 



Substituting this into Eq. (|A.12[) we get 

*W(ari <x 2 < ... < x N ) oc ^(-1) [P] ( f[ 



N r 



a<b 



q - -k{N + 1 - 2P a ) - lq--K(N + l-2P b )) - in 



x exp 



/V 



N 



iqJ2 x a + ^J2( N+1 - 2P *) X « 



/ N s r N N 

<* E(- 1 ) [P1 ( n t p " - p « + !] ) ex p "?E x - + \ + 1 - 2P ^ x - 

P ^a<b ' <- o=l a=l 



(B.13) 
(B.14) 



Here one can easily note that due to the presence of the product Il^<b [Pb — Pa . + 1] hi the summation over permutations 
only the trivial one, P a = a, gives non-zero contribution (if we permute any two numbers in the sequence 1,2, N 
then we can always find two numbers a < b, such that Pb = P a — 1). Thus 



vW(xi <x 2 < ... < x N ) 



oc exp 



N 



a=l 



a=l 



(B.15) 



Taking into account the relation, Eq. (|B.9[ ). we recover the function, Eq. (|B.l[) . which is symmetric with respect to its 
N arguments and therefore can be extended beyond the sector, Eq. (|A.2l) . for arbitrary particles positions. Finally, 
substituting the momenta, Eq. (|B.12l) . into the general expression for the energy spectrum, Eq. (|A.13p . we get 



a=l 



(B.16) 



Performing here simple summations (using Eqs. (|B.4p . (IB.5p ) one recovers Eq. (|B.6p . 



2. Eigenfunctions 



A generic eigenfunction of attractive bosons is characterized by momenta parameters {q a } (a = 1, 2, ...N) which 
may have imaginary parts. It is convenient to group these parameters into M (1 < M < N) "vector" momenta, 

q? = q a - \ k {n a + 1 - 2r) (B.17) 

where q a (a — 1,2,...,M) are the continuous (real) parameters, and the discrete imaginary components of each 
"vector" are labeled by an index r = 1, 2, n a . With the given total number of particles equal to N, the integers n a 
have to satisfy the constraint 

M 

n Q = N (B.18) 

a=l 

In other words, a generic eigenstate is characterized by the discrete number M of complex "vector" momenta, by the 
set of M integer parameters {ni, ri2, tim} = n (which are the numbers of imaginary components of each "vector") 
and by the set of M real continuous momenta {gi, q 2 , qui} = Q- 
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The general expression for the eigenfunctions is given in Eqs. (|A.9P - (|A.12p . To understand the structure of the 
determinant of the N x N matrix exp(iq a Xb) , which defines these wave functions, the N momenta q a , cq. (lB.17p . can 
be ordered as follows: 



{q a } = {q?} = {q{, <£, ... , q^; q\, cj, ... , ... ; gf , gf, ... , <} 



(B.19) 



By definition, 



det 



exp(iq a x c ) 



(c,d)=l,...,N 



^(-l) [p] exp\i^2q Pa x a 



(B.20) 



where the summation goes over the permutations of N momenta {q a }, Eq. (|B.19p . over TV particles {x\,X2, xjv}, and 
[P] denotes the parity of the permutation. For a given permutation P a particle number a is attributed a momentum 

component g"/"\ • The particles getting the momenta with the same a (having the same real part q a ) will be called 
belonging to a cluster fi a . For a given permutation P the particles belonging to the same cluster are numbered by 
the "internal" index r = 1, ...,n a . Thus, according to Eq. (|A.lip . 



/ N v r N 

^n'w = c'eh) |pi nh^-^)+ ! ' ts s n ( i «- 1 ')] exp *e 



a<b 



a(c) 
*r(c) Xc 



(B.21) 



where is the normalization constant to be defined later. Substituting here Eq. (|B.17p and taking derivatives we 

get 

N r 

'O*) = CM^-V^U (?«(-) "Y [na(a) + l-2r(o)])-(g a(6) -y[n a(6) +l-2r(6)])+i«sgn(x -a! 6 ) 



a<b 



x exp 



N N 

^2 ^(c)x c + ^ ^2 { n *(c) + 1 - 2r(c)^ x c 



(B. 



The pre-exponential product in the above equation contains two types of term: the pairs of points (a, b) which belong 
to different clusters (a(a) ^ ot(b)), and pairs of points which belong to the same cluster (a(a) = a(b)). In the last 
case, the product n a over the pairs of points which belong to a cluster fi Q reduces to 



n Q OC Yl [ r ( b ) ~ r ( a ) ~ S S n ( x a - Xb)]] 



(B.23) 



a<beO Q 



As for the ground state wave function Eq. (|B.14[) - (|B.15[) . one can easily note that due to the presence of this product 
in the summations over n Q ! "internal" (inside the cluster Q a ) permutations r(a) only one permutation gives non- 
zero contribution. To prove this statement, we note that the wave function ^/^''(x) is symmetric with respect to 
permutations of its N arguments {x a }; it is then sufficient to consider the case where the positions of the particles 
are ordered, x± < X2 < ■ ■ ■ < xn- In particular, the particles {x ak } (k — 1, 2, . . . , n a ) belonging to the same cluster 
Q a are also ordered x ai < x a2 < ■ ■ ■ < x ana ■ In this case 



n Q oc jj[r(Z)-r(*) + l] 



(B.24) 



k<l 



Now it is evident that the above product is non-zero only for the trivial permutation, r{k) = k (since if we permute 
any two numbers in the sequence 1, 2, ...,n a , we can always find two numbers k < I, such that r(l) = r(k) — 1). In 
this case 



n Q oc n[z-fc + i] 



(B.25) 



k<l 
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Including the values of all these "internal" products, Eq. (|B.25[) . into the redefined normalization constant C q , n , for 
the wave function, Eq. (IB.22p (with x\ < X2 < • • • < Xff), we obtain 

N p 

= CW^'(-1)11 1] (««00 - Y»«(a)) - (fc(6) " f »«*(*)) + i*(r(o) - r(6) - 1) 



x exp 



a<6 
a(a)^a(6) 

N 



N N 
i E Qa(c)%c + g E ( n «( c ) + 1 ~~ 2r ( c )) Xc 



(B.26) 



where the product goes now only over the pairs of particles belonging to different clusters, and the symbol means 
that the summation goes only over the permutations P in which the "internal" indices r(a) are ordered inside each 
cluster. Note that although the positions of particles belonging to the same cluster are ordered, the mutual positions 
of particles belonging to different clusters could be arbitrary, so that geometrically the clusters are free to "penetrate" 
each other. In other words, the name "cluster" does to assume geometrically compact particles positions. 

Now taking into account the symmetry of the wave function $qfn(x) with respect to the permutations of its 
arguments the expression in Eq. (|B.26l) can be easily continued beyond the the sector X\ < 22 < •■• < %n for the entire 
coordinate space R^. Using the relations 



^ ( n a + 1 - 2r(a))x a = ^(n a + I - 2k)x ak = \x ak 



(B.27) 



a£il a 



k=l 



k,l=l 



(where x ai < x a2 < ... < x an ), for the wave function "fq^n (x) with arbitrary particles positions we get the following 
sufficiently compact representation (cf. Eq. (|B.21[ l): 



a<b 
a(a)^a(b) 



-i(d Xa -d Xb ) +iKSgn(x a - x b ) 



exp 



M n a M n a 

: E 9« E *c - - E 



a=l cef2 Q 



a=l c,c'£f2 a 



Finally, substituting Eq. (|B.17|) - (|B.18j) into Eq. (|A.13|) . for the energy spectrum one easily obtains: 

M(q,n) = 27 E J2(ik) 2 = ^E - ^E( n « - n «) 



E 



a — 1 r—1 



2/3 



24/3 



(B.28) 
(B.29) 



a=l 



3. Orthonormality 

We define the overlap of two wave functions characterized by two sets of parameters, (M, n, q) and (M', n', q') as 



~.(M,M')t i\ 

QLn' (q.q) 



^x*K>)*W(x) 



Substituting here Eq. (IB.28|) we get 

O W) - C ) C"E E (- 1 ) H+|P1 / dN * 



(B.30) 



(B.31) 



p p> 



N 

n 

a<b 
a (a)^a' (b) 

N 

n 

a<b 
a(a)^a(b) 



i(d Xa - d Xb ) - iKSgn(x a - x b ) 



~i(dx a - d Xb ) + insgn(x a - x b ) 



M n 



exp 



exp 



M' n' a 

E s «-tE E \ x - 



a=l c£f2' 



M n„ 



a—1 c.c'gO' 



M n a 

1 E E Xc - 1 E E 1^ - »c 



q=i cen„ 



a—1 c,c'£fi Q 
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where {ft a } and {£l' a } denote the clusters of the permutations P and P' correspondingly. Integrating by parts we 
obtain 



P P' 



M' 



E ^ E Xc - 7 E E \ x c~ x c 



N 

n 

a (a)^a'(6) 



-i(9x„ - <9:r b ) - iKSgn(a; a - x b ) 



a=l c£iV a 
N 



a—1 c,c'Gfi' 



n 

a(a)y^a(6) 



-i(^x a - 9 Xb ) + iKSgn(x a - Xb) 



x exp 



M n a M n Q 

E E Xc - z E E 



a=l cen Q 



a=l c,c'eO Q 



(B.32) 



First, let us consider the case when the integer parameters of the two functions coincide, M = M', n = n', and 
for the moment let us suppose that all these integer parameters {n a } are different, 1 < ni < n-i < ... < um- Then, 
in the summations over the permutations in Eq. (IB.32[) . we find two types of terms: 

(A) the "diagonal" ones in which the two permutations coincide, P = P' ; 

(B) the "off-diagonal" ones in which the two permutations are different, P ^ P' . 
The contribution of the " diagonal" ones reeds 

/+oo r M n a M n a 

^xexp -^E^E^-iE E i*«-*c 
/ ~°° <- q=i cen„ q=i c.c'en„ 



N 

n 

a<b 
*(a)=£a(b) 



(9x a - 9x b ) + K Z 



exp 



c6Q 



M n a 



E E x ° - 7 E E 



3; c X c 



a—1 CjC'GOq 



(B.33) 



It is evident that all permutations a(a) in the above equation give the same contribution and therefore it is sufficient 
to consider only the contribution of the "trivial" permutation which is represented by line in Eq. (|B.19[) . The cluster 
ordering given by this permutation we denote by ao(a). For this particular configuration of clusters we can redefine 
the particles numbering, so that instead of a "plane" index a — 1,2, ...,N the particles would be counted by two 
indices (a, r): {x a } — > {xf} {a — 1,...,M) (r — l,...,n a ) indicating to which cluster a a given particle belongs and 
what is its "internal" cluster number r. Due to the symmetry of the integrated expression in Eq. (|B.33P with respect 
to the permutations of the particles inside the clusters, we can introduce the "internal" particles ordering for every 
cluster: xf < x% < ■■■ < x" .In this way, using the relation, Eq. (IB.27p . we get 



)(M,M) {A) 



Nl 



ni\n%\...nM 



M 

n 



-\-oo 



dx°t 



dxl 



x exp 



M n a 



a — 1 r—1 

M n a np f 



M n a 



E^E< + ?EEK +1 - 2r ) x r 

M n a 

E««E< + ? EE( n « + 1 - 2r >r 



nnn -(<w-^r+« 

a</3 r—1 r' — 1 L 



exp 



a—1 r—1 



M n a 



(B.34) 



where the factor N\/nx\...nj\#\ is the total number of permutations of M clusters over N particles. Taking the 
derivatives and reorganizing the terms we obtain 



M n a ™p 



{q a - y»o) - (qp - +in(r-r' - 1) 



« m)<a W) = cmcm; m n n n 

^a</3 r=l r'=l 

| Z' + OO Z' + OO Z' + OO 1 

xJJW d:E ? / dx%.... d< o e' ( «°^» )£ "°i I ° +fiE °=i (n ° +1 - 2r)l ° (B.35) 

a=l I J x ? ^; o -i J 
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Simple integrations over x™ yields (cf. Eqs. (IB.8l) - (|B.10l) ) 



« M)(A, (q,q') = i^r^dinn 



M n a nf, 



M 

n 



(n a \) 2 K 



a</3 r=l r'=l 

(2n)S(q a - q' a 



(q a - y n ") ~ {ip - y"/ 3 ) +*«(*"- r' - 1) 



(B.36) 



Now let us prove that the "off-diagonal" terms of Eq. (IB.32p . in which the permutations P and P' are different, 
give no contribution. Here we can also chose one of the permutations, say the permutation P, to be the "trivial" 
one represented by line in Eq. (|B.19[) with the cluster ordering denoted by ao(a). Given the symmetry of the wave 
functions it will be sufficient to consider the contribution of the sector x\ < X2 < ■■■ < %n- According to Eq. (|B.32l) . 
we get 



. r M n a M n a 

« M)<B W) oc E'(-i) [p,] / exp -ij^q'* E *a-zE E I*. 

P' Jx 1 <...<x N I a=1 Qgn , a=la,bEQ' 



Xb 



N 

n 

a<b 
a (a)^a (b) 



-i{d Xa - d Xb ) +iKSgn{x a - x b ) 



N 

n 

a<b 
a' {a)^a ! (b) 



-i(d Xa -d Xb ) - iKSgn(x a - x b ) 



x exp 



M n a M n a 

' E Qa E Xa + 9 E E + 1 _ 2r {°)) x a 



Q=l a6f2° 



(B.37) 



Here the symbols {fi^} denote the clusters of the trivial permutation ao(a). Since P' ^ P, some of the clusters ft' a 
must be different from f2° . As an illustration, let us consider a particular case of N — 10, with three clusters ri\ = 5 
(denoted by the symbol "O") ; n 2 = 2 (denoted by the symbol " x") and = 3 (denoted by the symbol "A"): 



particle number a 


1 


2 


3 


4 


5 


G 


7 


8 


9 


10 


permutation ao(a) 

















X 


X 


A 


A 


A 


permutation a' (a) 











A 





X 


X 





A 


A 



Here in the permutation a' (a) the particle a = 4 belong to the cluster a = 3 (and not to the cluster a = 1 as in 
the permutation ao(a)), and the particle a = 8 belong to the cluster a = 1 (and not to the cluster a = 3 as in the 
permutation ao(a)). Now let us look carefully at the structure of the products in Eq. (|B.37[) . Unlike the first product, 
which contains no "internal" products among particles belonging to the cluster Q°, the second product does. Besides, 
the signs of the differential operators {d Xa — d Xb ) in the second product is opposite to the "normal" ones in the first 
product (cf. Eqs. l|B.23[) - (|B.25jl ). It is these two factors (the presence of the "internal" products and the "wrong" 
signs of the differential operators) which makes the "off-diagonal" contributions, Eq. (|B.37l) . to be zero. Indeed, in the 
above example, the second product contains the term 



K 5 = 



-i(d Xi - dx 5 ) + in 



exp 



3 n a 

z E^ E 

a— 1 a£Q% 



\ E E ( n « 

a— 1 a£EO° 



1 - 2r(a))x a 



(B.38) 



(we remind that the particles in the clusters f2° are ordered, and in particular x$ < £5) 
get 



Taking the derivatives, we 



n' 



4.5 



ki + |0*i + 1 - 2r(4)) - iqi - |(m + 1 - 2r(5)) 



[r(4) - r(5) + l] = 



(B.39) 



since in the first cluster r(a) = a. 
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One can easily understand that the above example reflect the general situation. Since all the cluster sizes n a are 
supposed to be different, whatever the permutation a'(a) is, we can always find a cluster f2° such that some of its 
particles belong to the same cluster number a in the permutation a 1 (a) while the others do not. Then one has to 
consider the contribution of the product of two neighboring number points 



k,k+l 



-i(d Xk -d Xk+1 ) +in 



exp 



M 



E 



M 



a— 1 a£Q,° 



(B.40) 



where in the permutation a! (a) the particle k belong to the cluster number a and the particle (fc + 1) belong to some 
other cluster. Taking the derivatives one gets 



K, k+ i * [r(fc)-r(fc + l) + l] = 
as r(a) is the "internal" particle number in the cluster fi°, where r(k + 1) 







(B.41) 



r(Jfc) + 1 (cf. Eas. (lR^ - (lB~25l) '). 
Thus, the only non-zero contribution to the overlap, Eq. (|B.30[) . of two wave function ^7*2 ( x ) an< ^ ^o!* 1 ( x ) (having 
the same number of clusters M and characterized by the same set of the integer parameters 1 < n\ < n2 < ... < tim) 
comes from the "diagonal" terms, Eq. (|B.36P : 



M 



\\2 K n a 



M n a 



nnn 

. V</ir=lr'=l 



in (r 



1) 



M 

n 



(2n)S(q a - q' a ) 



(B.42) 



The situation when there are clusters which have the same numbers of particles n a is somewhat more complicated. 



Let us consider the overlap between two wave function ^gt^(x) and \?q,n^(x) (which, as before have the same M 



,( M ), 



and n) such that in the set of M integers n\,n2, ■■■,tim there are two n a 's which are equal, say n ai 



(where 



ot\ 7^ ot-i). In the eigenstate (q',n) these two clusters have the center of mass momenta q' ai and q' a2 , and in the the 
eigenstate (q, n) they have the momenta q ai and q a2 correspondingly. According to the above discussion, the non-zero 
contributions in the summation over the cluster permutations a(a) and a' (a) in Eq. (|B.32[) appears only if the clusters 
{f2 Q } of the permutation a(a) totally coincide with the clusters of the permutation a'(a). In the case when 

all n a are different this is possible only if the permutation a(a) coincides with the permutation a' (a). In contrast to 
that, in the case when we have n ai = n a21 there are two non-zero options. The first one, as before, is given by the 
"diagonal" terms with a(a) = a' (a) (so that the clusters {tt a } and {£l' a } are just the same), and this contribution 
is proportional to 5(q ai — 9a 1 )^(9a 2 — Qa 2 )- ^ ne second ("off-diagonal") contribution is given by such permutation 
a' (a) in which the cluster Q,' (of the permutation a' (a)) coincide with the cluster tt a2 (of the permutation a(a)) and 
the cluster Q! a2 (of the permutation a' (a)) coincide with the cluster Q ai (of the permutation a(a)) while the rest of 
the clusters of these two permutations are the same, Cl' a = Cl a (a ^ 0:1,0:2). Correspondingly, this last contribution 
is proportional to S(q ai — q' a2 )S(q a2 — q' ai ) (— l) n< *i. In fact this situation with two equivalent contributions is the 

consequence of the symmetry of the wave function ^^^(x): the permutation of two momenta q ai and q a2 belonging 
to the clusters which have the same numbers of particles, n ai = n a2 produces the factor (— l)" a i. This is evident 
from the general expression for the wave function, eq. (|A.9p . where the permutation of any two momenta q ai and q a2 
belonging to the clusters which have the same numbers of particles corresponds to the permutation of n columns of 
the matrix exp{iq a Xb) ■ Therefore considering the clusters with equal numbers of particles as equivalent and restricting 
analysis to the sectors q ai < q a2 ; q' < q' a2 we find that the second contribution, S(q ai ~q' a2 ) 5{q a2 ~Qa 1 ) ^ s identically 
equal to zero, thus returning to the above result Eq. (IB.42p . 

A generic eigenstate (q, n) with M clusters could be specified in terms of the following set of parameters: 



(q> n ) = {(gi,mi),...,(q Sl ,mi);(g Sl+ i,m 2 ),...,(g Sl+S2 ,m 2 ); 



(Qsi 



+ ...+Sfc-l + l 



Wfc),..., (q Sl+ ... +Sk ,m k )} (B.43) 



where si + s 2 + ... + Sk = M and k integers {mi} (1 < k < M) are all supposed to be different: 

1 < mi < m 2 < ••• < rrik 
Here the integer parameter k denotes the number of different cluster types. For a given k 

M k 



zli Ua = X/ Simi = N 



(B.44) 



(B.45) 
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Due to the symmetry with respect to the momenta permutations inside the subsets of equal n's it is sufficient to 
consider the wave functions in the sectors 

qi < q 2 < ... < q Sl ; (B.46) 

?si+i < q Sl +2 < ••• < q Sl +s 2 ; 



<?»l+...+a*-l+l < q Sl +...+s k _ 1 +2 < ••• < <7ai+...+a*-l+a* 

In this representation we again recover the above result Eq. (|R42|) 

Finally, let us consider the overlap of two eigenstates described by two different sets of integer parameters, n' ^ n. 
In fact this situation is quite simple because if the clusters of the two states are different from each other, it means 
that in the summation over the pairs of permutations P and P' in Eq. (|B.32l) there exist no two permutations for 
which these two sets of clusters {£l a } and {Fl' a } would coincide. Which, according to the above analysis, means that 
this expression is equal to zero. Note that the condition M' ^ M automatically implies that n' ^ n. 

Thus we have proved that 



Cn M 'W) = \C^\ 2 8(M,M')(l[6( 



[n a ,n c 



M 



H(2n)5(q a -q' a ; 



<m n 



M n a 



n n n i ( ? ° y n °) ~ ~ + iK ^ - r ' ~ ^ 

. a</3 r=l r' = l 



(B.47) 



where the integer parameters {n a } and {n' a } are assumed to have the generic structure represented in Eqs. ([B.43jl - 
(|B.45|) . and the momenta {q a } and {q' a } of the clusters with equal numbers of particles are restricted in the sectors, 
Eq. (|B.46p . According to Eq. (|B.47l) . the orthonormality condition defines the normalization constant 



r M 

n 



(n a \ 



TIq, fx 



M n a 



nnnr — — — — 

a <fs r=i r'=i (q a - ^n a ) - (q fj - fn f3 ) + in (r - r' - 1) 



(B.4 



In other words, the wave functions, Eq. (IB.28p . form the orthonormal set. Although, at present we are not able to prove 
that this set is complete, the suggestion of the completeness (which assumes that there are exist no other eigenstates 
besides those described above) looks quite natural. 



4. Propagator 

The time dependent solution ^(x, t) of the imaginary-time Schrodinger equation 

N N 

/9fit*(x;t) = + -Kj2^(x a -x b )^( X ;t) (B.49) 

a—1 a^b 

with the initial condition 

tf(x;0)=Ilf =1 <5«) (B.50) 
can be represented in terms of the linear combination of the eigenfunctions ^^/(x), Eq. (IB. 28ft : 

*(*,<) = J2 E' /' Pq ^n(x)*^ r (0) exp[-£ M (q,n)t] (B.51) 

M=l n ^ 

where the energy spectrum Em{<\, n) is given by Eq. (|B.29p . The summations over n Q are performed here in terms of 
the parameters {s.;,mi}, Eqs. (|B.43P - (|B.45l) : 

E' s E E E 5 (E s " M ) *(E s ^> (b.52) 

n fc=l si...Sfe=l l<mi...<m fc M=l ' M=l ' 
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where 8(n,l) is the Kronecker symbol, and for simplicity (due to the presence of these Kronecker symbols) the 
summations over m, and are extended to infinity. The symbol J Pq in Eq. (IB.51|) denotes the integration over M 
momenta q a in the sectors, Eq. (|B.46|) . 

The replica partition function Z{N, L) of the original directed polymer problem is obtained via a particular choice 
of the final-point coordinates, 



v / r' 

Z(N,L) = *(0;L) = ]T ]T Vq |*^(0)| 2 exp[-E M (q,n) L] 

1\ IT 1 — « 

According to Eq. (|B~28l) . for M > 2, 



(B.53) 



M=l n 



*^(o) = ^)^'(-i) [pi n 



a<b 
cz(a)^cx(b) 



to(a) - 17 [n-a(a) + 1 - 2r(a)] ) - ( q a ( b ) - — [n a (jt,) + 1 - 2r 



"" m!n 2 !...n M ! 



M n„ 



n n n ~ y^) _ (*> - y^) + ^ - r o 



a</3 r=l i-'=l 

Substituting here the value of the normalization constant, Eq. (|B.48[) . we get 



(B.54) 



|«(o)| 2 = 



NIk 



N 



M 

n 



n;=i nXi (<?* - f n Q ) - & - % np ) + iK ( r - o 



[Ua=i Kn <*) c<0 UZi UXi I (to - t™«) - (to ~ T»vO + i« (r - - 1) 
This expression can be essentially simplified. Shifting the product over r' in the denominator by 1 we obtain 



(B.55) 



|«(0) 



N\k 



N 



Ua=l KU <* ) £<P Ur=l (to - f »«) - (?/? - f »/») + »« (r - ^ - 1) 



M 

n 



n"=r 



2 ' ''CK 



) - (?/s - f n 0) +*«(»■- 1) 



(B.56) 



Redefining the product parameter r in the denominator, r — > n a + 1 — r, and changing the obtained expression (under 
the modulus square) by its complex conjugate we get 



1*^(0) 



NIk 



N 



Jla=l Kn «) «<0 Ur=l (q a ~ fn a ) - (qp - fnp) 
Shifting now the product over r in the numerator by 1 we finally obtain 



M 

n 



K=i (to - ^n a ) - {qp ~ f np) +m{r- 1) 



(B.57) 



1*^(0) 



N\k n 



M 

n 



IL=l ra M a <p | to -9/3 ~ f ("a + n /3)| 



For M = 1, according to Eqs. (lB~TT) and (|B~TT|) . 



1*^(0)1 



k n N\ 
kN 



(B.58) 



(B.59) 



Since the function /(q, n) = l^q^n (0)| 2 exp[— Em(<1, n) L] in Eq. (|B.53l) is symmetric with respect to permutations 
of all its M pairs of arguments (q a , n a ) the integrations over M momenta q a can be extended beyond the sector defined 
in Eq. (|B.46p for the whole space Rm- As a consequence, there is no need to distinguish equal and different n Q 's any 
more, and instead of Eq. (|B.52l) . we can sum over M integer parameters n a with the only constrain, Eq. (|B.18|) (note 
that this kind of simplifications holds only for the specific "zero final-point" object ^(0; t), Eq. (|B.53p , and not for the 
general propagator ^(x;i), Eq. (|B.51[) containing N arbitrary coordinates x\, ...,xn)- Thus, instead of Eq. (|B.53[) we 
get 



Z(N,L) 



N r M r+00 , CO -, / Af v 

E i n / f E l«(o)| 2 e 



-E M (q,n)L 



(B.60) 
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Substituting here Eqs. (|B.29p . (|B.58|) and (|B.59I) we get the following sufficiently compact representation for the replica 
partition function: 



Z(N.L) — Nl k 



N 



-(- oo 



dq 



exp 



, 2wkN 

N ^ r AI oo f+c , 

e w\ n e 



Ml 

M=2 l q=1 n„=l 



dq a 



x exp 



4E n ^ Q + ^E( n « - n «) 

M a=l f q=1 



M 



oo 

n 2 L 



2irKn c 

M 



M 



a</3 



|?« - 9/3 - f ("a + rtp)\ 



(B.61) 



The first term in the above expression is the contribution of the ground state (M = 1), while the next terms (M > 2) 
are the contributions of the rest of the energy spectrum. After simple algebra the above replica partition function 
can be represented as follows: 



z(n,l) =e 



-fSNLf 



Z(N,X) 



(B.62) 



where / = ^(3 4 u 2 — j^ln(f3 3 u), and 



Z(N.L) = Nl 



+oo 



dq 



exp 



2itkN 

M oo „+oo 



2(3 q 24/3 



+ e jjt n n 

M=2 L ct=l n a = l 



dq 



x exp 



L 



— oo 

M k 2 L M il 



2jTKn c 

M 



5 E""^ n 



a</9 1 9a -9/3 - Tf(»l Q + na)| 



24/3 ^ 

a = l 



(B.63) 



Appendix C 

Fredholm determinant with the Airy kernel and the Tracy- Widom distribution 



In this Appendix the original derivation of Tracy and Widom [14[ will be repeated in simple terms to demonstarate 
that the function F2 (s) defined as the Fredholm determinant with the Airy kernel can be expressed in terms of the 
solution of the Panleve II differential equation, namely 



F 2 (s) = det[l - K A ] = exp 



dt{t - s)q 2 (t) 



where Ka is the integral operator defined on semi-infinite interval [s, 00) with the Airy kernel, 

Ai(t 1 )Ai'(t 2 )-Ai'(t 1 )Ai(t 2 ) 

K A (t 1 ,t 2 ) = 

ii — t 2 

and the function q(t) is the solution of the Panleve II equation, 

q" =tq + 2q 3 

with the boundary condition, q{t — ) +00) ~ Ai(i). 



Let us introduce a new function R(t) such that 

F 2 (s) = exp 



dtR{t) 



(C.l) 

(C.2) 
(C.3) 



(C.4) 
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or, according to the definition, Eq. (|C.ll) . 



R( s ) = — ln|det(l 



Here the logarithm of the determinant can be expressed in terms of the trace: 
In 



det(l-K A ) = -J2~ Tr ^A 

71=1 

= dt i dt 2 ... dt n K A (t 1 ,t 2 )K A (t 2 ,t 3 )...K A (t n ,t 1 ) 

n=l U •* s ^ s ^ s 



(C.5) 



(C.6) 



Taking derivative of this expression we gets 
R(s) = - 



dt(l-K A ) (s,t)K A (t,s) 

00 pOO pOO pOO 

= -K A (s,s)-J2 / dt l / *2- / dt n - 1 K A (8,ti)K A (t 1 ,t2)...K A (t n - 1 ,s) 

Substituting here the integral representation of the Airy kernel, Eq. (|C.2l) . 

K A (h,t 2 ) = / dz Ai(h + z) Ai(i 2 +2) 
Jo 

after some efforts in simple algebra one gets 



(C.7) 



R(s) 



00 p 00 



dti / d* 2 Ai(ti) (1 - if a) (*i,t 2 ) Ai(t 2 



(C.8) 



(C.9) 



Taking the derivative of this expression and applying some more efforts in slightly more complicated algebra, we 
obtain 



ds 



R(s) = -q 2 (s) 



where 



According to Eq. ([CT0)) , 



Let us introduce two more functions 



q(s) = / dt(l-K A ) \s,t) Ai(t) 



R(s) = / dtq\t) 



(CIO) 



(C.ll) 



(C.12) 



poo poo 

v(s) = J dt x J dt 2 Ai(t 1 )(l-K A y 1 (t 1> t 2 ) Ai'(t 2 ) 

poo 

p{s) = / dt(l-K A )~\8,t) Ai'(i) 



(C.13) 
(C.14) 

Taking derivatives of the above three functions q(s), v(s) and p(s), Eqs. (|C.ll[) . (jC.13l) and (|C.14[) . after somewhat 
painfull algebra one finds the following three relations: 



q' = p - Rq 

p = s q — p R — 2qv 

v = — p q 



(C.15) 
(C.16) 
(C.17) 
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Taking derivative of the combination (R 2 — 2v) and using Eqs. (|C.10|) and (|C.17|) . we get 

^-(R 2 -2v) = 2q(p - Rq) 
On the other hand, multiplying Ea. (IC.15p by 2q we find 



(C.18) 



(C.19) 



Comparing Eqs. (|C.18P and (|C.19|) and taking into account that the value of all the above functions at s — > oo is zero, 
we obtain the following relation 

R 2 - 2v = q 2 (C.20) 
Finally, taking the derivatiove of Eq. (|C.15|) and using Eas. (|C.10|l . (|C.15[) . (|C.16[) and (|C.20|) we easily find 

q" = 2q 3 + sq (C.21) 



which is the special case of the Panleve II differential equation [36|, |37|. Thus, substituting Eq. (|C.12[) into Eq. (|C.4p 
we obtain Eq. (|C.ll) . 

The function -Fifs) gives the probability that a random quantity t described by a probability distribution functions 
Prw(t) has the value bigger than a given parameter s: 



dtP TW (t) 



Taking the derivative of this relation and subtituting here the result, Eq. (|C.ip , we find 



Ptw(s) = exp 



dt(t - s)q 2 (t) 



dtq 2 (t) 



(C.22) 



(C.23) 



In the limit s — > oo the function q(s), according to its definition, eq. ljC.lip . must go to zero, and in this case 
Eq. (|C.2ip turns into the Airy function equation, q" = sq. Thus 



q(s oo) ~ Ai(s) ~ exp 



s 3/2 



(C.24) 



It can be proved |38j that in the opposite limit, s — >• — oo, the asymptotic form of the solution of the Panlevee equation 
(|C.2ip (which has the right tail Airy function limit, Eq. (|C.24IO is 



q(s — > — oo) 



1 

'2 S 



(C.25) 



Substituting the above two asymptotics into Eq. (|C.23l) . we can estimate the asymptotic behaviour for the right and 
the left tails of the TW probability distribution function: 



Ptw( s +oo) ~ exp 
Ptw( s — oo) ~ exp 



i s 3/ 2 



12' 



(C.26) 
(C.27) 
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